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' Abstract 

>Y 

■ In this work, we investigate the dynamics of a non-local model describing spontaneous 

cell polarisation. It consists in a drift-diffusion equation set in the half-space, with the 
coupling involving the trace value on the boundary. We characterize the following be- 
haviours in the one-dimensional case: solutions are global if the mass is below the critical 
mass and they blow-up in finite time above the critical mass. The higher-dimensional case 
is also discussed. The results are reminiscent of the classical Keller-Segel system in double 
Oh ■ the dimension. In addition, in the one-dimensional case we prove quantitative convergence 

results using relative entropy techniques. This work is complemented with a more realistic 
model that takes into account dynamical exchange of molecular content at the boundary. 
In the one-dimensional case we prove that blow-up is prevented. Furthermore, density 
converges towards a non trivial stationary configuration. 
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1 Introduction 

> 

Cell polarisation refers generically to a process that enables a cell to switch from a spheri- 
cally symmetric shape to a state with a prefered axis. Such a phenomenon is an essential 
■^j- \ step for many biological processes and is involved for instance in cell migration, division, or 

morphogenesis. While the precise biochemical basis of polarisation can vary greatly, in its 
early stages polarisation is always characterised by an inhomogeneous distribution of specific 
molecular markers. Cell polarisation can be driven by an external asymmetric signal as in the 
example of chemotaxis, where a chemical gradient imposes the direction of migration of cells 
[2]. Another example is given by mating yeast, for which the external signal is a pheromone 
^ \ gradient, which causes the cell to grow an elongation known as a shmoo in the direction of 

the pheromone source [2]. However observations show that some cellular systems, such as 
mating yeast, can also polarise spontaneously in absence of external gradients [34J. These 
two distinct polarisation processes, driven or spontaneous, are necessary for cells to fulfil dif- 
ferent biological functions. However, so far the conditions under which a cell can polarise 
spontaneously or only in response to an external asymmetric forcing are not well understood. 
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The molecular basis of cell polarisation has been much discussed in the biological literature 
over the past decade, and is likely to involve several processes. It is now widely recognised 
that the cell cytoskeleton plays a crucial role in cell polarisation. The cell cytoskeleton 
is a network of long semiflexible filaments made up of protein subunits (mainly actin or 
microtubules). These filaments act as roads along which motor proteins are able to perform a 
biased ballistic motion and carry various molecules, in a process which consumes the chemical 
energy of adenosine triphosphate ATP. It is observed that the efficiency of formation of polar 
caps in yeast, indicating polarisation, is reduced when actin transport is disrupted, and that 
the polar caps formed are unstable [34, 35, 22J. In the case of neurons, it has been shown 
that the polarisation of the growth cone is suppressed when microtubules are depolymerised 
[8]. To account for these observations, it is generally argued that the cytoskeleton filaments 
mediate an effective positive feedback in the dynamics of polarisation markers |34J. This arises 
from the molecular markers not only diffusing in the cell cytoplasm, but also being actively 
transported by molecular motors along cytoskeleton filaments, the dynamic organisation of 
which is regulated by the markers themselves. 

From the physical point of view, achieving an inhomogeneous distribution of diffusing 
molecules without an external asymmetric field as in the case of spontaneous polarisation 
requires either an interaction between the molecules or a driving force that maintains the 
system out of equilibrium. In the case of the cell cytoskeleton, it is well known that the 
hydrolysis of ATP acts as a sustained energy input which drives the system out of equilib- 
rium, and one can therefore hypothesizes on general grounds that spontaneous polarisation 
in cells stems from non equilibrium processes. Cell polarisation has been the subject of a 
few theoretical studies in recent years. Many models rely on reaction-diffusion systems in 
which polarisation emerges as a type of Turing instability [231 [251 ES] and some (e.g. [2TJI [54] ) 
include cytoskeleton proteins as a regulatory factor. However, the full dynamics of markers 
is generally not considered. 

In this article, following the work of [TBJ, we study a class of models for spontaneous cell 
polarisation. These models couple the evolution of molecular markers with the dynamics of 
the cytoskeleton. Namely the markers are assumed to diffuse in the cytoplasm and to be 
actively transported along the cytoskeleton. The density of molecular markers is denoted by 
n(t,x). The advection field is denoted by u(t,x). This field is obtained through a coupling 
with the boundary value of n(i, x). 

The cell is figured by the half-space % = M^ -1 x (0, +oo). We denote the space vari- 
able x = (y,z). The time evolution of the molecular markers follows an advection-diffusion 
equation: 

d t n(t,x) = An(t,x) - V • (n(t,x)u(t,x)) , t>0, xeU. (1.1) 
1.1 The one-dimensional case 

We first analyse two different models set on the half-line (0, +oo). In the simplified version, 
the advection field is given by u(t, z) = — n(t, 0). Active transport arises at uniform speed, 
the speed being given by the value of the density at z = 0. 

1.1.1 The simplified model 

The model writes as follows. 

d t n(t,z) = d zz n(t,z) +n(t,Q)d z n(t,z) , t > , z e (0, +oo) , (1.2) 
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together with the zero-flux boundary condition at z = 0: 

d z n(t,0) +n(t,0) 2 = 0. (1.3) 
We have formally conservation of molecular content: 



M= n (z)dz= / n(t,z)dz. 

Jz>0 Jz>0 

Solutions of (|1.2p may become unbounded in finite time (so-called blow-up). This occurs 
if the mass M is above the critical mass: M > 1. In the case M < 1, the solution converges 
to 0. In the critical case M = 1 there exists a family of stationary states parametrized by 
the first moment. The solution converges to the stationary state corresponding to the first 
moment of the initial condition J z>0 zno(z) dz. 

Theorem 1.1 (Global existence and asymptotic behaviour in the sub-critical and critical 
cases: M < 1). Assume that the initial data tiq satisfies both tiq £ + z) dz) and 

Iz>o n o( z )0-°S n o( z ))+ dz < +oo. Assume in addition that M < 1, then there exists a global 
weak solution (in the sense of Definition \2.1\) that satisfies the following estimates for all 
T > 0, 



sup / n(t, z)(log n(t, z)) + dz < +oo , 
te(o,T) Jz>o 

/ / n(i, z) {d z \ogn{t,z)) 2 dzdt < +oo . 
Jo Jz>0 



In the sub-critical case M < 1 the solution strongly converges in L 1 towards the self-similar 
profile G given by \2.26\) in the following sense: 



lim 

t— >>+oo 



n(t,z) 1 G 



VTT2i \VTT2i 



= 0. 

L 1 



In the critical case M = 1, assuming in addition that the second moment is finite f z>0 z 2 no(z) dz < 
+oo, the solution strongly converges in L 1 towards a stationary state aexp(— az), where 
a.- 1 = f z>Q zn (z)dz. 



Theorem 1.2 (Blow-up of weak solutions: M > 1). Assume M > 1. Any weak solution with 
non-increasing initial data no blows-up in finite time. 



In the present biological context, blow-up of solutions is interpreted as polarisation of the 
cell. Indeed there is a strong instability driving the system towards an inhomogeneous state. 

In Section [3l we present analogous blow-up results in the case of a finite interval z G (0, L) 
or finite range of action. 

Remark 1. Such a critical mass phenomenon (global existence versus blow-up) has been widely 
studied for the Keller-Segel system (also known as the Smoluchowski-Poisson system) in two 
dimensions of space [30] . The equation (jl.2p represents in some sense a caricatural version 
of the classical Keller-Segel system in the half-line (0, +oo). Note that there exist other ways 
to mimick the two dimensional case in one dimension |12|, [T3] . 



3 



Remark 2. There is a strong connection between the equation under interest here (jl.2p and 
the one-dimensional Stefan problem. The later writes indeed [19) : 

( d t u(t, z) = d zz u(t, z), t > , z e (-oo, s(t)) , 

\ lim^-oo 8 z u(t, z)=0, u(t, s(t)) = , d z u(t, s(t)) = -s'(t) . 

The temperature is initially non-negative: u(0,z) = u${z) > 0. By performing the following 
change of variables: (p(t, z) = —u(t,s(t) — z), we get an equation that is linked to (|2.ip by 
n(t,z) = d z (p(t,z). This connection provides some insights concerning the possible continu- 
ation of solutions after blow-up [19]. This question has raised a lot of interest in the past 
recent years [201 E21 EH HI] ■ It is postulated in [19] that the one-dimensional Stefan problem 
is generically non continuable after the blow-up time. 

1.1.2 The model with dynamical exchange of markers at the boundary 

The boundary condition (|1.3p turns out to be unrealistic from a biophysical viewpoint. This 
claim is emphasized by the possible occurence of blow-up in finite time. On the way towards a 
more realistic model, we distinguish between cytoplasmic content n(t, z) and the concentration 
of trapped molecule on the boundary at z = 0: fi(t). Then the exchange of molecules at the 
boundary is described by very simple kinetics: 

^/x(t) = n(t,0)-7M(i). 
The transport speed is modified accordingly: u(t, z) = —fj,(t). The model writes: 

d t n(t, z) = d zz n(t, z) + fj,(t)d z n(t, z) , t > , z G (0, +oo) 
d z n(t,0) + n(t)n{t,0) = 

The flux condition on the boundary ensures the conservation of molecular content. De- 
noting m{t) = j z>0 n(t, z) dz the partial mass of cytoplasmic markers, we have: 

M = no + m = fi(t) + m(t) . 

Since the transport speed is bounded, fj,(t) < M, we clearly have global existence of so- 
lutions for any mass M > 0. We can precise the asymptotic behaviour in the super-critical 
case M > 1. This is the purpose of the following Theorem. 

Theorem 1.3. Assume that the initial datano satisfies bothno E L 1 ((l+z) dz) and f z>0 no(z)(log no(z))+ dz 
+oo. Assume the mass is super- critical M > 1. The partial mass m(t) converges to 1 and 
the density n(t,z) strongly converges in L 1 towards the exponential profile (M — l)e~( M ~ 

1.2 The higher-dimensional case 

In the higher dimensional case N > 2 we only partially analyse simplified models such as 
(jl.2p where the transport speed is directly computed from the trace value n(t, y, 0). Equation 
(jl.ip is complemented with the zero-flux boundary condition: 

d z n(t,y,0) - n(t,y,0)u(t,y,0) ■ e z = 0, i/eK^ 1 . (1.4) 
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We have formally conservation of the molecular content: 



M = I uq(x) dx 
Jh 




n(t, x) dx . 



Following |18| we make the distinction between two possible choices for the advection 
speed u. In the transversal case, the field u is normal to the boundary: 



This corresponds to a particular orientation of the cytoskeleton, modelling the microtubules. 
Indeed microtubules are very rigid filaments whose bending length is larger than the typical 
size of yeast cells. 

In the potential case, the field u derives from a harmonic potential. The source term of 
the potential is located on the boundary: 



This corresponds to another orientation of the cytoskeleton, modelling the actin network. 
Indeed the actin networks is a diffusive network where orientations are mixed up. In dimension 
N = 1, observe that the two choices (|1.5p and (jl.6p coincide. 

In dimension N > 2, we state global existence for small initial data. The criteria are iden- 
tical for the two possible choices of the advection field (|l,5p or (jl.6j) . This is a consequence 
of the two common features: both fields are divergence free and possess the same normal 
component at the boundary. 

Theorem 1.4 (Global existence in dimension N > 2). Assume that the advection field sat- 
isfies the two following conditions: V • u > and u(t,y,0) ■ e z = n(t,y,0). Assume that the 
initial data uq satisfies both no E L l ((l + \x\ 2 ) dx) and ||71o[|ljv *s smaller than some constant 
cn depending only on the dimension N . Then there exists a global weak solution to and 



Notice that both conditions V • u > and u(t, y, 0) • e 2 = n(t, y, 0) are fulfilled in (jl.5p 
and (fL6j) . 

Theorem 1.5 (Blow-up in dimension N > 2). Assume that n(t,x) is a strong solution to 
which verifies: 

• d z n(t, x) < for all x £ % and t > when the advective field is given by il.5\) . 

• d z n(t,x) < and for all x £ % and t > 0, the matrix A(t,x) = x (g> d z V y log n(t, x) 
satisfies A T + A > (in the matrix sense) when the advective field is given by lil.b]) . 

Assume in addition that the second momentum is initially small enough: there exists a con- 

N+l 

stant Cm depending only on the dimension such that J x£ ^\x\ no(x) dx < C^M . Then 
the maximal time of existence of the solution is finite. 



u(t,y,z) 



n(t,y,0)e z . 



(1.5) 




(1.6) 



{HP- 
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Open questions We end this introductory Section with some open questions that we are 
not able to resolve, (i) Obtain a rate for the convergence in relative entropy in Theorem 11.11 
for the cases M = 1 and M < 1. (ii) Prove blow-up for the systems (|l.ip - (|1.6p with large 
initial data without any monotonicity assumption on the density n(t,x). 

The outline of the paper is as follows. In Section [21 we analyse with full details the 
one-dimensional case. In section [3] we study some variants of blow-up criteria in the one- 
dimensional case. In Section [U we study a model with flux of markers at the boundary in 
the one-dimensional case. In Section [5j we analyse the higher dimensional case. 

Results in the one-dimensional case have been announced in the note 1111. 



2 The boundary Keller-Segel (BKS) equation in dimension 

N = 1 

In this Section we study the following equation, 

f d t n(t, z) = d zz n(t, z) + n(t, 0)d z n(t, z) , t > , z G (0, +oo) , 
[ d z n(t,0) +n(t,0) 2 = 0, 

and we prove Theorems 11.11 and 11.21 More precisely, in Sections 12.11 we prove the existence of 
a global weak solution for M < 1. Then in Section [2. 31 we prove the blow up character in the 
case M > 1. 

We begin with a proper definition of weak solutions, adapted to our context. 
Definition 2.1. We say that n(t, z) is a weak solution of (|2.ip on (0, T) if it satsifies: 

n G L°°(0,T;LV(M + )), d z n G L 1 ((0, T) x M + ) , (2.2) 
and n(t, z) is a solution of (|2.ip in the sense of distributions in Z>'(R_|_) 

Since the flux (d z n(t,z) + n(t,0)n(t, z)) belongs to L 1 ((0, T) x R+), the solution is well- 
defined in the distributional sense under assumption (|2.2p . In fact we can write n(t, 0) dt = 

~ Jo L>o d Mt, z) d z dt - 

Weak solutions in the sense of Definition 12.11 are mass-preserving: 



M= n (z)dz= / n(t,z)dz. 

Jz>0 Jz>0 

The proof closely follows the arguments of the next Lemma which is concerned with moment 
growth. 



Lemma 2.2 (Moment growth). Assume n(t,z) is a weak solution of Assume in addi- 

tion that zuq G L 1 (M + ). Then the following identity holds true: 

[ zn(T,z)dz= [ zn (z)dz+ [ (l-f n(t, z) dz) n(t, 0) dt . (2.3) 

Jz>0 Jz>0 JO V J z>0 J 
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Proof. Consider the approximation function x( z ) which verifies x( z ) = 1 if < 2 < 1, 
x( z ) = if z > 2, which is smooth and non-negative everywhere. Define the family of 
functions (<f e ) e by (f e ( z ) = z x( £z )- We recall the weak formulation: 



n(T, z)(p £ (z) dz = I n (z)ip £ ( z )dz 

[ (d z n(t, z) + n(t, 0)n(t, z)) <p' e (z) dz dt . 

Jz>0 

The function (f £ ( z ) converges monotically to z as e — > 0, hence from the monotone convergence 
theorem, we deduce that zn(T,z) £ L . 

The function (f' E (z) = x( ez ) + ez x' ( £z ) 1S bounded in L°° uniformly in e and it converges to 
1 a.e. Since n(-,0)n G L 1 ((0, T) x R + ) and d z n G L 1 ((0, T) xR + ), from Lebesgue's dominated 
convergence theorem, it follows that 



£-5-0 



fT p i-T i- 

lim / / ip' £ (z)n(t,0)n(t,z)dzdt = n{t,0)n(t, z) dz dt , 

" '0 Jz>0 JO Jz>0 

lim / I ip' £ (z)d z n(t,z)dzdt = [ [ d z n(t, z) dz dt = - [ n(t,0)dt. 

" '0 Jz>0 JO Jz>0 Jo 



£-S>0 



□ 



2.1 Global existence for sub-critical mass M < 1 
2.1.1 A priori estimates 

Our next result is concerned with the derivation of a priori bounds for solutions to (|2.ip in 
the classical sense. 

Proposition 2.3 (Main a priori estimate). Let n be a classical solution to \2.1\) . If M < 1, 

then the following estimate holds true for some 5 > and for all t G (0,T): 

n(t, <z)(log n(t, z)) + dz + 5 / / n(s, z)(d z logn(s, z)) 2 dzds 
z>o Jo Jz>0 

<[ n (z) (log n (z)) + dz+ f zn (z) dz + C(T) . (2.4) 

Jz>0 J z>0 

Proof. We first derive the following trace-type inequality. 

2 

\2 ' ' 



n(t,Qy = / d z n(t,z)dz 



z>0 



< (/ n(t,z)dz)( n(t,z)(d z \ogn(t,z)y dz) . (2.5) 

I z>0 J \Jz>0 



We compute the evolution of the entropy 

— I n(t, z) log n(t, z) dz = ( dtn(t, z) log n(t, z) dz 
di J 2>o J z >o 



z>0 Jz>0 

(d z n(t, z) + n(t, 0)n(t, z)) ^J^lA dz 
z>o n{t,z) 

n(t, z) (d z log n(t, z)) 2 dz + n(t, 0) 2 . (2.6) 

2>0 
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The two contributions are competing. We estimate the balance using inequality (12.51) . 



_d 



n(t, z) log n(t, z) dz < (M — 1) / n(t, z) (d z log n(t, z)) 2 dz . 

z>0 Jz>0 



On the contrary to the classical two-dimensional Keller-Segel equation, the dissipation of en- 
tropy gives directly the sharp criterion on the mass. There is no need to seek a free energy 
as in [7] (and references therein). To control the negative part of the entropy, we use the 
following Lemma adapted from [10] . 

Lemma 2.4. For any f G Li.(K+, (1 + z)dz), if f /log / < +oo ; then / log / is in L 1 (]R + ) 
and for all a > 0, the following inequality holds true: 

[ /(z)(log /(*))+ dz< / f (z) (log f(z) + az) dz + — . (2.7) 

Jz>0 Jz>0 ae 

Proof. Let / = /l/<i and m = f z>Q f(z)dz. We build up the relative entropy between / 
and ae~ az . 

[ J(z) (log7(*) + az) dz = [ log (J^L) ae— dz + mloga . 

Using Jensen's inequality, we deduce that 

lo jm W „ d2 



>0 ae az \ae 



>( / -I^Lae~ az dz\loK ( ! -I^-ae- az dz 



z>o ae az ) \J Z>0 ae az 



Therefore, 



Using 



= m log m . 

I f(z) log f(z) dz + a I zf(z) dz > m log (am) > 



f(z)(\ogf(z)) + dz= [ f(z)logf(z)dz- [ f(z)logf(z)dz, 

z>0 Jz>0 Jz>0 

this completes the proof of Lemma 12.41 □ 
Let us now estimate the first moment. Recalling (|2.3|) . we deduce that 

zn(t,z)dz < / zuq(z) dz + / n(s,0)ds 
z>0 Jz>0 Jo 



< [ zn (z)dz + ^- + 5' f n(s,0) 2 ds, 
Jz>o 4 5 Jo 

< / zn (z)dz + — 
Jz>o 4c" 



2>0 

+5' I I n(s,z) (d z log n(s,z)) 2 dzds. (2.8) 

Jo J z>0 
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Combining f|2T6j) . (j2T7|) and (|2T8j) with a = 1 we obtain that 



n(i, z) (log n(t, z)) + dz + (1 — M — 5') / / n(s, z) (d z log n(s, z)) 2 dz ds 

n (z)logn (z)dz + / zn (z)dzH 

2 >o 7 2 >o e 4<5' 



< 



Since M < 1 we can choose 5' > such that (|2.4p holds. □ 
2.1.2 Regularization procedure 

To prove existence of weak solutions in the sense of Definition 12.11 we perform a classical 
regularization procedure. We carefully choose our function spaces in order to end up with 
minimal assumptions on the initial data. We introduce 



(*) = / Mz)n £ (t,z)dz, 

Jz>0 



where 4> £ is an approximation to the identity. We have formally a £ (t) — > n(t, 0) as e — > 0. 
We consider the following regularized problem 

f d t n £ (t, z) = d zz n s (t, z) + a £ (t)d z n'(t, z) , 

{ 2.9) 
I d z n £ (t,0) + a £ (t)n £ (t,0) = 0. 



Our aim is to extend the main a priori estimate (|2.4p to the regularized problem (|2.9p . We 
check that 



a £ {t) = - <f> e (z) d z n £ (t,y)dydz< \d z n £ (t, z)\ dz . 

Thus the following inequality replaces (|2.5|) : 



a £ (t)n £ (t,0) < Ml / n £ (t,z)(d z log n £ (t,z)) z dz 

\Jz>0 

On the other hand the moment growth estimate only relies on the diffusion contribution. We 
have accordingly, 

f zn £ (t,z)dz < [ zn (z)dz + ^- (2.10) 
Jz>o Jz>o 4c» 

+5' I [ n £ (s,z) (d z log n £ (s, z)) 2 dzds. 

JO Jz>0 

It is then straightforward to justify (|2.4p for the regularized solution n £ in the line of Propo- 
sition [2]3j There exists 5 > such that 

n £ (t,z) (log n £ (t,z)) + dz + 5 f f n £ (s, z)(d z logn e (s, z)) 2 ds dz 

z>0 JO J z>0 

< n (z) (log n (z)) + dz+ I zn (z) dz + C (T) . (2.11) 

Jz>0 Jz>0 
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2.1.3 Time compactness 



Passing to the limit as e — > 0, the main difficulty lies in the nonlinear term a £ (t)d z n £ (t, z). 
We need some compactness to proceed further. It is provided by the Aubin-Simon Lemma, 
see [261 EI]. 



Lemma 2.5 (Aubin-Simon). Let X C B C Y be Banach spaces such that the embedding 
X C B is compact. Assume that the set of functions J- satisfies: J- is bounded in L 2 (0,T; X) 
and dtf is uniformly bounded in L 2 (0,T;Y). Then T is relatively compact in L 2 (0, T;B). 



The natural choice for spaces in our context would be X = W rl - 1 (R + ) and B = C°(M+) 
(up to the decay problem at infinity). However, due to the possible apparition of jumps, the 
embedding X C B is not compact. Using the entropy estimate (|2.1ip we are able to modify 
the space X in order to make the embedding X C B compact. The crucial point is to obtain 
an equi-continuity condition weaker than any Holder condition, in the spirit of pQ Theorem 
8.36]. 

Lemma 2.6. Assume T is a set of non-negative bounded functions in the following sense: 
there exists a constant A > such that for all f G T 

sup / f(t, 2) (log f(t, z)) + dz < A, [ T [ f(t,z)(d z logf(t,z)) 2 dzdt< A. 

te(0,T) Jz>0 JO Jz>0 

Then there exists a continuous function r/ : M + — > R + and a constant A' depending on A such 
that r](0) = and for all function f G T we have 

2 

I dt<A'. (2.12) 



r ( \m 

/ sup 

JO \x^y 



y)-f(t,x)\ 



i](y - x) 

Proof. First for x < y we have that 



2 



\f(t,y)-f(t,x)\ 2 < U \d z f(t,z)\dz 

< (j" f(t,z)dz^j^J(t,z)(d z logf(t,z)) 2 dz^ . (2.13) 
We use the Jensen's inequality for x < y: 

' 1 V f(t, z) dz) log (— r f(t, z) dz) < — r f(t, z)Qog f(t, z))+ dz 



y-xj x J \y~xj x J + y-x 

A 

< 



y-x 



We can invert this inequality to get: 

— ! V fit, z)dz<<S> (—) , (2.14) 

y-xJx \y-xj 

where $ : [0, +00) — > [1, +00) is the reciprocal bijection of x(logx)+. We define 77(2) = 
z<&iz~ l A). Clearly rjiz) — > as z — > since $ is sublinear. Combining (|2.13|) and (|2.14p . we 
deduce the estimate (j2.12j) . □ 
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We denote C ' 71 the space of functions having modulus of continuity controlled by r\: 
C « = lg€C°: sup ^ - g( f }l < +ool . 

{ x+ y my-x) J 

The injection C 0,v C C° is compact on bounded intervals [T]. The behaviour of functions 
outside bounded intervals in our context is controlled by the following estimate which is a 
consequence of (|2.13|) as y — > +00: 



l/(*,^)r <-U zf(t,z)dz \ U f(t,z)(d,]agf(t,z)y dz \ . (2.15) 

The last requirement in the Aubin-Simon Lemma consists in getting very weak estimate 
for the time derivative dtn £ . We can write 

d t n £ (t,z)+d z f(t,z) = 0, 

where j £ (t,z) = d z n £ (t,z) + a £ (t)n £ (t, z) is uniformly bounded in L 2 (0,T; L 1 due to 

(12. lip and the following inequalities: 

IKHi2 (0iT) < l|dX||£2 (0jT ;£i(R + )) ^ M ^ J^n £ (t,z)(d z \ogn £ {t,z)) 2 dtdz. 

Hence <9 t n £ is uniformly bounded in L 2 (0,T; (W rl '°°(R + )) / ) . 

We introduce some useful functional spaces, endowed with their corresponding norms: 

X = {geC°"(R + ): z l ' 2 g{z) e L™{^ + )} , 
B = C (M + )nL°°(R + ), 
Y = (W 1 ' QO (R+))' . 

It is straightforward to check that X is compactly embedded in B. 

Combining the above estimates ()2.10I - [2~TT1 - I2.12ll2.15p we obtain that n £ is bounded in 
L 2 (0,T; X) uniformly with respect to e. The Aubin-Simon Lemma ensures that, up to ex- 
tracting a subsequence, n £ converges strongly in L 2 (0,T; B) towards some n. Prom uniform 
convergence of n £ , we deduce that a £ (t) — > n(t,0) strongly in L 2 (0, T). Hence we can pass to 
the limit in the nonlinear term a £ (t)n £ (t, z) in the weak formulation. 

To conclude we verify that the a priori estimates given in Proposition 12.31 are valid after 
passing to the limit e — > 0. From the strong convergence in L 2 (0,T; B) we deduce that, up 
to extracting a subsequence that we do not relabel, 

lim / n £ (t, z) (logn e (t, z)) + dz = I n(t, z) (log n(t, z)) + dz , a.e. t 6(0, T) . 

£ ^°Jz>0 ' Jz>0 

On the other hand, we use the convex character of the functional (see [7] and the references 
therein) j z>Q f(z) (d z log f{z)) 2 dz = 4 j z>Q (d zy /f(z)j dz. We have finally, 

liminf / / n £ (s, z) (d z log n £ (s, z)) 2 dz ds > / / n(s, z) {d z log n(s, z)) 2 dz ds . 
Jo J z >o Jo Jz>0 

So the a priori estimate fj2.4[) is valid a.e. t G (0,T). 



11 



2.2 Long-time behaviour for the critical and the subcritical cases 

In this Section, we investigate long-time behaviour of solutions in the case M < 1 using 
entropy methods. We distinguish between the critical case (no need to rescale) and the 
sub-critical case (self-similar diffusive scaling). 

We stress out that the method for proving global existence in the critical case M = 1 
strongly relies on the entropy estimate (|2.16p . This is why we analyse the global existence 
and the long time behaviour all in all. 

2.2.1 The critical case: global existence and asymptotic convergence 

The main inequality we have used so far in order to prove global existence is (|2.5p . Equality oc- 
curs if logn(i, z) is linear w.r.t. z: there exists a(t) > such that n(t, z) = Ma[t) exp(— a(t)z). 
In fact the boundary condition (|2.ip implies M = 1. On the other hand the stationary states 
to equation (|2.1|) are precisely the one-parameter family: 

h a (z) = aexp (—az) , a > . 

This motivates to introduce the relative entropy: 

H(t) = f ^f^log J h a {z)dz = [ n(t, z) log n(t, z) dz + aJ(t) -logo. 

Jz>0 ha{Z) V J J z> 

Recalling (|2.3|) . we notice that the first momentum of density is conserved in the case 
M = 1: J(t) = J(0). This prescribes the value for a provided we can pass to the limit t — > oo: 
a -1 = J(0). We also recall the formal computation giving the time evolution of the relative 
entropy (|2.6p : 

— H(t) = -/ n(t,z)(d z log n{t,z)) 2 dz + n(t,0) 2 (2.16) 

at Jz>0 

= - n(t,z) (d z logn{t,z) +n(t,0)) 2 dz < 0. 

Jz>0 

The Jensen's inequality yields H(t) > 0, so we have < H(t) < H(0). We deduce from 
Lemma 12.41 that the quantity J z>0 n(t, 2) (log n(t, z)) + dz is uniformly bounded by some con- 
stant denoted by Cq: 

I n(t, z)(log n(t, z)) + dz < Cq , a.e. t G (0, +00) . 

Jz>0 

The method of proving convergence in relative entropy towards h a is as follows. We first 
gain a priori estimates which enable to pass to the limit after extraction as in Section 12.1.31 
For this we update the estimates in Section [2. II with the key information that the entropy H 
is uniformly bounded. The identification of the limit requires more information concerning 
the behaviour of the density at infinity. We use the fact that the first momentum drives the 
evolution of the second one. Finally we conclude that the entropy converges to along some 
subsequence. Since it is non-increasing, it converges to globally. 
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A-priori bound We cannot follow the strategy developped in Section \2. II since we crucially 
used M < 1. We need to gain some control on the dissipation (|2.16p which is the competition 
of two opposite contributions that are nearly equally balanced. For that purpose we introduce 

1 /2 

the function A : R + -4 R + such that A(0) = and A' (it) = (logu)^ . It is non-decreasing, 
convex and superlinear. Thus there exists A £ R such that A(u) 2 > 2Cqv? for all it > A. 
Adapting (|2.5p to our context we get 

A(n(i,0)) 2 



< 
< 
< 

From (I2TT61) and (l2T7j) . we deduce that 

d f u M < + \a <r J ° if n M)^ A ' 

— / n(t, z) log n(t, z) dz < < A(n(t,0)) 2 . U ^ U n\2 -r U n\ >s A 
dtj z>0 — Co +n(t,0y < -n{t,0) z if n(t,0)>A. 

We introduce the set E = {t : n(t,0) > ^4}. We have obtained the estimate 

/ n(t,0) 2 dt< / n {z)logn (z)dz, (2.18) 

JE Jz>0 

thus n(t, 0) cannot be too large (in I? sense). 

We deduce from (|2.16p and f|2.18j> that f f z>0 n(s, z) (d z logn(s, z)) 2 dzds is bounded for 
all t G (0, T). The previous statements prove that Proposition 12.31 remains valid in the case 
M = 1. Next, the existence proof, in the case M = 1, is similar to the case M < 1 and we do 
not repeat it here. 



d z (A(n(t,z))) dz 



2>0 



>0 



A'(n(t, z))ra(i, z)5 z (log n(i, z)) dz 



2>0 



n(t,z)|A'(n(t,z))| 2 dz J ( / n(t, z) (d z log n(t, z)) 2 dz 



2>0 



Z>0 

C / n(t,z) (9 2 log n(t,z)) 2 dz . 

Jz>0 



n(t,z) (log n(t,z)) + dz) [ / n(t, z) (d z log n(t, z)) 2 dz 



Z>0 



(2.17) 



Passing to the limit Let A r be any integer. We translate the solution in time: we define 
ujsr(s, x) = n(N + s,x). The function H(t) is non-increasing and bounded below by zero. 
Therefore the entropy dissipation (|2.16j) converges to zero in an averaged sense. The estimate 



f + ( [ n(t,z)(d z logn(t,z)) 2 dz-n(t,0)A dt = H(N)-U(N+l) 
JN \Jz>0 J 



N~>oo 



->0. 



reads 

»l 



un(s, z) (d z log u N (s, z)) 2 dz -u N (s, 0) 2 ds > . (2.19) 



7V->oo 



'0 \Jz>0 

We deduce from ()2.18p that ujv(s,0) is bounded in L 2 (0, 1) uniformly w.r.t. N. Hence both 
terms are bounded in (|2.19|) . This enables to pass to the limit as in Section 12.1.31 Up to 
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extracting a subsequence (labelled with N') there exists Uoo such that u^' u oo strongly in 
L 2 (0,l;5): 

[ \\u N ,(s)- Uoo (s)\\ 2 B ds^0. (2.20) 
Jo 

We can pass to the limit in each term of the averaged dissipation: 

= liminf / I / un'(s, z) (d z logit/v'( s > z)) 2 dz — un'(s, 0) 2 ) ds 



o \Jz>0 



N'-+oo 
-l 

^ j. „. /„ n^2 



> / / ^oo(s,z) (d z logUoo(s,z)) dz-iioo(s,0) ds>0. 

J \Jz>0 J 

We have used the L 2 (0, T; L°°(M. + )) strong convergence (|2.20p to pass to the limit in the 
nonlinear term un'(s,0) 2 , and also the convexity of the functional J z>0 f(z) (d z log f(z)) 2 dz 
(see Section l2,1.3p . 

Identification of the limit We deduce that Moo satisfies almost everywhere 

u QO (s, z) = f3(s) exp(— a(s)z) , a(s), (3(s) > . 

To determine a(s) and f3(s) we shall use the conservations of mass and first momentum. 
Since the first momentum is uniformly bounded, we have that M = lim f z>0 v>N'(t, z) dz = 
Iz>o u oo(t, z) dz. This yields a(s) = /3(s). 

We have proved so far that we can always extract a subsequence such that un'(s,z) 
approaches u OQ (s,z) in L 2 (0,T; B). We explain below why it is delicate to derive a{s) = 
a = J(0) _1 without any better control of the density n(t, z) as z — > +oo. Suppose we have 
a(s) = a and the convergence un(s, z) — > Uoo(z) is uniform. We would have on the one hand, 

a _1 = liminf / zuN'(t, z) dz > / zuoo(z) dz = (a) -1 , (2-21) 

Jz>0 Jz>0 

and on the other hand, 

0<limH(t)= / u co (z) log Uoo(z) dz + 1 — log a = log a — log a . 

Jz>0 

We would deduce a > a which is the same as (|2.21[) . 

In the case J z>Q z 2 uq(z) dz < +oo, let us examinate the evolution of the second momen- 
tum. We simply have 

Id/* 

-— z 2 n(t,z)dz = M -n(t,0)3(t) = l-n(t,0)a- 1 . (2.22) 
2 dt J z>0 

The idea is to pass to the pointwise limit n(t, 0) — > a. If a > a, the right-hand side of (|2.22p 
becomes asymptotically 1 — aa -1 < which leads to a contradiction. 
Let introduce the notation I(t) = f z>Q (z 2 /2)n(t,z)dz. We have 

1{N + 1) - I(JV) = / (1 - n(t, 0)a _1 ) dt = / (l - u N (s, 0)a _1 ) ds . 

J N Jo 
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Since I is a non-negative quantity, we clearly have limsupI(A r + 1) — i(N) > 0. Furthermore 
we have lim sup I(JV + 1) — I(iV) < 0. To see this, assume on the contrary that lim sup I(N + 
1) — I(A r ) = 5 > 0. We can extract a converging subsequence. Keeping the same notations as 
above, we have limI(iV' + 1) — I(iV') = 5. We can pass to the limit similarly (up to further 
extracting) in the following average quantities: 



a 1 = liminf / / zujft(t, z) dzds 
Jo Jz>0 



> f f zu oa (s,z)dzds = [ (a(s)y 1 ds, (2.23) 

Jo Jz>0 JO 

5 = lim/ (l-u N >(s,0)a~ l ) ds = 1 - of 1 f a(s)ds. (2.24) 
Jo Jo 

Inequality (|2.23p yields a(s) ds > a by Jensen's inequality. This is in contradiction with 

(MP - 

We conclude that lim sup I(N + 1) — I(JV) = 0. We extract a converging subsequence, such 
that limI(iV' + 1) - I(N') = 0. Hence we obtain (pT23D and (I2T241 with 5 = 0. The equality 
case in Jensen's inequality yields a(s) = a. 

Asymptotic convergence (without any extraction) We have proved that there exists 
a subsequence such that Ujv 7 converges towards Uqq — ha in L 2 (0,T;B). We cannot pass 
to the limit pointwise in time from L 2 convergence. However there exists a sequence of 
times sn> £ (0,1) such that \\un'{sn') — ^oo||b - > 0. This includes uniform convergence 
and uniform decay at infinity. We can pass to the limit in the entropy term and we obtain 
H(un'(sn')) H(uoo) = 0. It means H[n(N' +sjy>)] — > 0. Using the non-increasing property 
of the entropy we have H[n(t)\ — > as t — > oo (without extracting any subsequence). 

Finally we recall the Csiszar-Kullback inequality [151 12"3] . For any non- negative functions 
/, g G L^lL).) such that J x>0 f(x) dx = f x>0 d(x) dx = 1, the following inequality holds true, 

||/- 9 ||| 1(Rt) <4/ i>o/ (x),o g (M) dx . 
This yields \\n(t) — h a \\ L i — > 0. 

2.2.2 Self-similar decay in the sub-critical case 

In the sub-critical case M < 1 the density n(t, z) is expected to decay with a self-similar 
diffusion scaling [7]. To catch this asymptotic behaviour we rescale the density accordingly: 

n(t,z) = n(logy / l + 2t, — F === ) . 

The new density u(t, y) satisfies: 

d T u(r, y) = d yy u(T, y) + d y (yu(r, y)) + n(r, 0)d y u(T, y) , (2.25) 

together with a no-flux boundary condition: d y u(r, 0) + u(t, 0) 2 = 0. The additionnal left- 
sided drift contributes to confine the mass in the new frame (r,y). The unique stationary 
equilibrium in this new setting can be computed explicitely: 

G a (y) = aexp (-ay - y 2 /2) , (2.26) 
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where a is uniquely defined by the condition J y>0 G a (y) dy = M. This rewrites P(a) = M, 
P being an increasing function defined as follows: 

P{a) = L ^ h - &) ^ ' { lim-To! P(a) = 1 • 
We re-define the relative entropy and the first momentum in the rescaled frame: 



J(r) = / yu(r,y)dy. 
Jv>0 



>y>0 

We also introduce a Lyapunov functional for equation ()2.25|) : 

L(T) = H(T) + 2(1 — M) (J(T) " a(1 " M))2 • 
Note that it is a non-negative quantity by Jensen's inequality. 
Lemma 2.7. The Lyapunov functional L is non-increasing: 

JL L (t) = -D(t)<0. 

The dissipation reads as follows 

D(t) = / u(r,y) (d y log u(t, y) +y + u(r,0)) 2 dy 

Jy>0 

+ (T^M)(i JW ) 2 ' (2 ' 27) 



Proof. We compute the evolution of the entropy as previously: 
-^H(r) = J d T u(r,y) ^log(u(r,y)) + ay + y^j dy 

—, — V + a + y) d y 



y>0 

u ( T > y) {dy lo g u ( T > y) + yf dy + u ( T > °) 2 - u ( T > °) J ( r ) 

y>0 

+an(r,0) - aJ(r) - au(r,0)M 

U ( T > y) lo g u ( T > y) + y + o)) 2 dy 

y>0 

+(M - l)u(r, 0) 2 + u(r, O)J(r) + a(l - M)n(r, 0) — aJ(r) . (2.28) 
Moreover, the time evolution of the first momentum reads in the rescaled frame: 

-^J(r) = (l-M)n(r,0)-J(r). 
dr 
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As compared to (|2.32[) the additional contribution is due to the rescaling drift. We can 
eliminate u(r, 0) from (|2.28|) in the two following steps: 

u(T,0)J(r)+a(l-Jtf)u(r,0)-aJ(r) = J^Jf^ (r) + + a^-3(r) 

2 



and 



d (J(t) — a(l — M)) 
~ dr" 2(1 - M) 

2J(r) d JM , J(r) 2 , . 

+ (T3M)d7 J(r) + (r3M)' (2 - 29) 



4-3(t)) 2 = (M - l)«(r, 0) 2 + J*(rL jL j( r ) + . (2.30) 



(1 — M) \ dr / v / v > / (1 — M) dr (1 — M) 

Combining (|2T2"5]) - (|2T2"§|) - ([230]) the proof of Lemma Ellis complete. □ 

To prove convergence of u(t, ■) towards G a we develop the same strategy as in Section r2.2.1l 
for the critical case M = 1. The main argument (apart from passing to the limit) consists in 
identifying the possible configurations for which the dissipation D vanishes. In fact this 
occurs if and only if both terms in (|2.27p are zero. This means that Joo(r) = (1 — M)uoo(r, 0) 
on the one hand, and on the other hand, 

dy\ogUoo{T,y) + y + iioo(r,0) = 0. 

We obtain that Uqo = G Q , where G a is given by ()2.26p . To pass to the limit as in Section 
I2.2.1l we need to gain some good control of J y>Q u(t, y) (d y logn(r, y)) 2 dy from the dissipation 
term D. The situation here is simpler than in Section [2.2.11 since the mass is sub-critical. The 
argument goes as follows 

U ( T , U) {dy log u(t, y)+y + n(r, 0)) 2 dy 

y>o 

u(t, y) {d y log n(r, y)) 2 dy + (M — 2)«(r, 0) 2 + 2u(r, 0)J(r) 



y>o 

+ [ y 2 u(T,y)dy-2M 

Jy>0 



> (m + -^-2Ju(t,0) 2 -2M, 

where we have used inequality (j2.5|) . The quantity M + M~ l — 2 is positive since M < 1. 
Hence, recalling Proposition 12.31 we can prove directly that u(-,0) belongs to I? locally in 
time (this was the purpose of (|2.17p - (|2.18|) ). 

Finally, we obtain that L converges to zero as r — > +oo. So u(r, •) converges towards G a 
in entropy sense. 

2.3 Blow-up of solutions for super-critical mass 

To prove that solutions blow-up in finite time when mass is super-critical M > 1 and no is 
non-increasing, we show that the first momentum of n(t, z) cannot remain positive for all 
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time. This technique was first used by Nagai [27], then by many authors in various contexts 
(see [51 El QU HSl E] for instance). 

The assumption that no is a non-increasing function guarantees that n(t, •) is also non- 
increasing for any time t > due to the maximum principle. In fact the derivative v(t, z) = 
d z n(t, z) satisfies a parabolic type equation without any source term, it is initially non-positive, 
and it is non-positive on the boundary due to (jl.3p . 

Therefore — d z n(t, z)/n(t, 0) is a probability density at any time t > 0. We deduce from 
the Jensen's inequality the following interpolation estimate: 

-d z n{tz) d \\ f j-Orttz) ^ 
z>o n (^0) ) J z>0 n(t,0) 

It rewrites in a more convenient way as follows, 

M 2 <2n(t,0) zn(t,z)dz. (2.31) 

Jz>0 

We denote the first momentum J(t) = j z>Q zn(t, z) dz. We plug ()2.3ip into the evolution 
of the moment (|2.3|) : 



J(t) = 3(0) + (1 - M) [ 

Jo 



t 

n(s, 0) ds 



, . (1 - M)M 2 f* 1 , . 

< J( o) + L_^y o _ ds . (2 . 32 ) 

We introduce the auxiliary function K(i) = J(0) + (1 — M)M 2 J(s) -1 ds. It is positive and 
it satisfies the following differential inequality: 

d K ^ _ (1 - M)M 2 1 < (1 - M)M 2 1 



dt w 2 J(t) - 2 K(t) ' 

hence, 

-^K(t) 2 < (l-M)M 2 . 

We obtain a contradiction: the maximal time of existence T* is necessarily finite when M > 1. 
On the other hand, following [23], it can be proved that the modulus of integrability has to 
become singular at T*: 



lim sup / (n(t, z) - K)+ dz) > . 



v iG(0,T*) Jz>0 

Otherwise a truncation method enables to prove local existence by replacing n with (n — K )+ 
for K sufficiently large. 

Remark 3. It is natural to perform the Laplace transform on the equation (|2.ip C z (n(t, z)) = 
n(t,Q = J z>0 n(t, z) exp(— Qz) dz. Then the occurence of blow-up is clear after transforma- 
tion. We refer the reader to [9] where the Fourier transformation has been applied successfully 
to analysing a one-dimensional caricature of the two-dimensional Keller-Segel equation. 
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3 Variants of blow-up criteria 



In this section we determine necessary conditions for blow-up to occur for a fast decaying 
interaction potential (Section 13. ip and for a finite interval (Section 13. 2p . 

3.1 Finite range of action 

In this part we consider the following system: 

d t n(t,z) = d zz n{t,z)-d z {n{t,z)dz<j>(t,z)) , t > , z e (0, +oo) , (3.1) 
with zero-flux at z = and the attractive potential is given by 

- d zz (j)(t, z) + a 2 (j)(t, z) = , -d z <t>(t,0) =n(t,0). (3.2) 
We introduce the exponential moment of the solution: 

Ja(t) = / exp(az)n(i, z) dz . 
Jz>0 

Proposition 3.1. Assume M > 1 and the exponential moment is small in the sense of 
criterion \3. 3\) below. Assume in addition that exp(— az)no(z) is a non-increasing function. 
Then the solution to \3. 1\) - \3. 2\) with initial data n(0,z) = uq{z) blows-up infinite time. 

Proof. The attractive field is given by d z (p(t,z) = — exp(— az)n(t, 0). Similarly to the proof 
of Theorem 11.21 we compute the time derivative of J a (t): 

— J a (t) = a 2 3 a (t) + an(t, 0)(1 - M) . 

We check that the function u(t, z) = exp(— az)n(t, z) is decreasing w.r.t. z for all time t > 0. 
For this purpose we write the equation for v(t,z) = d z u(t,z). This reads as follows 

dtu(t, z) = d zz u(t, z) + 2ad z u(t, z) + a 2 u(t, z) + exp(— az)n(t, 0)d z u(t, z) , 
dtv(t,z) = d zz v(t,z) + 2ad z v(t,z) + a 2 v(t,z) + exp(—az)n(t,0)d z v(t,z) 
—a exp(— az)n(t, 0)v(t, z) . 

Since the boundary condition reads v(t, 0) = —an(t, 0) — n(t, 0) 2 < and the above parabolic 
equation preserves non-positivity we deduce that v(t,z) < if v{0, z) < 0. 
We can adapt the inequality (|2.31|) to the function u(t, z) and we obtain 

M 4 < f [ exp(az)n(t, z) dz] ([ u(t,z)dz 

\Jz>0 J \Jz>0 

< J a (t) 2 n(t,0) 2 ( [ z ~ d f {t n Z) dz"" 

< j a (t)Mt,o) 2 [ f ^±LL^ )z^A dz 

Jz>o\ 2a 2 J u{t,0) 

< —J a (t) 2 n(t,0) / (exp(az) — exp(— az)) n(t, z) dz 
a Jz>o 



1 / M 2 

< -J a (t) 2 n(t,0)[J a (t)-—- 
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Finally, when M > 1 we obtain that: 



d T . . 2t . , a 2 (l — M)M 4 

-Mt) < a 2 3 a (t) + 



<(M-1). (3.3) 



Notice that J a (0) > M by definition. We get an obstruction to global existence if the following 
condition holds true, 

Jq(0) 4 ( _ JVf_ 2 
V Ja(0) 2 

□ 

3.2 Finite interval 

In this part we consider the equation (|1.2j) on a finite interval (0, L) for some L > 0, namely, 

d t n(t,z) = d zz n(t,z) + (n(t,0) - n(t,L))d z n(t,z) , t > , x € (0, L) , (3.4) 

together with ra(t = 0, z) = no(z) > and zero-flux boundary conditions at both sides of the 
interval. 

Equilibrium configurations are given by the family of functions: 

h{z) = aexp(— (a — 0)z) , j3 = aexp(— (a — f3)L) . (3-5) 

There are two possibilities, either a = j3 and h is constant, or a //3andM = f^h[z)dz = l. 
Observe that given a > there exists a unique /3 satisfying (|3,5p . If aL < 1 then f3 > a (h is 
increasing), whereas if aL > 1 then /3 < a (h is decreasing). 

Proposition 3.2. Assume M > 1 and i/ie /zrst moment is small: 4J(0) < LM. Assume in 
addition that n${z) is a non-increasing function. Then the solution to {3.$ with initial data 
n(0,z) = no(z) blows-up infinite time. 

Proof. We proceed again as in the proof of theorem (jl.2p . From Jensen's inequality, it follows 
that: 

-d z n(t,z) \ f L 2 -d z n(t,z) 

Z—, — — ; rrdz < / Z —, — — ; dz , 



l n(t, 0) — n(t, L) J ~ J n(t, 0) — n(t, L) 
hence, using that n(t, 0) > n(t, L) for any time t > 0, we deduce that 

(M — Ln(t, L)) 2 < (n(t, 0) — n(t, L)) ^2 f zn(t, z) dz - L 2 n(t, L)) , (3.6) 



o 

and the inequality remains true when n(t, 0) = n(t, L) and n(t, ■) is constant. Therefore, the 
first momentum J(t) = J Q L zn(t, z) dz satisfies: 

-^J(i) = (l-M)(n(t,0)-n(t,L)) 

< (1-M)^- Ln ^ 2 



2J(t) - L 2 n(t,L) 
M 2 - 2MLn(t, L) 



< (1-M) 

~ y 1 2J(t) 
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On the other hand, from (|3.6p again, it follows that 23 (t) > L 2 n(t, L) and we deduce that 

d „ , M(l-M) / 4J(t) 
-3(t) < — K , ' M 



dt w _ 2J(t) V L 
and the result follows by contradiction as in Section 12,31 □ 

4 The model with with dynamical exchange of markers at the 
boundary: prevention of blow-up and asymptotic behaviour 

In Section 12.31 we proved that finite blow-up occurs in the basic model (|2.ip when mass is 
super-critical M > 1. On the other hand the model which was originally proposed in [18] is 
the following: 

{d t n(t, z) = d zz n(t, z) + fi(t)d z n(t, z) , t > , z G (0, +oo) 
£ tt i(t) = n(t,0)-(i(t), 
together with the flux condition at the boundary: 

0X*,o) + M*M*,o) = ^(t)- (4.1) 

The quantity [i represents the concentration of markers which are sticked to the bound- 
ary and thus create the attracting drift. The dynamics of \i is driven by simple attach- 
ment/detachment kinetics. The mass of molecular markers is shared between the free parti- 
cles n(t, z) and the particles on the boundary fi(t). The boundary condition (|4.ip guarantees 
conservation of the total mass: 



n(t, z) dz + n(t) = M . (4.2) 



z>0 



Prom (|4.2p . we easily deduce that finite time blow-up cannot occur since the drift /i(t) is 
bounded by M. We denote by m(t) the mass of free particles: 



The conservation of mass reads 



m(t) = / n(t, z) dz . 

Jz>0 



We re-define the relative entropy as follows: 

where the asymptotic profile h is given by: 

h(z) = v exp {—vz) , u = M — 1 . 

When mass is super-critical M > 1, we shall prove that the density of free markers n(t,z) 
converges in relative entropy towards h, whereas the concentration of markers sticked at the 
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boundary fj,(t) converges to v. This is achieved using a suitable Lyapunov functional as in 
Sections 12.2.11 and 12.2.21 We introduce accordingly 

L(t) = m(t)H(t) + ~ (ji(t) - vf + fi(t) log + m(t) log m{t) . 

The rest of this Section is devoted to the proof of the following Lemma. 



Lemma 4.1. The Lyapunov functional L is non-increasing: 

^L(t) = -D(t)<0. 

TTie dissipation reads as follows 

D(t) = ^n(t,z)(a 2 logn(t,z) + ^M) 2 dz + m(t) - ' 

+ (n(t, 0) - /i(i)) log ( ! ^) + (MW " • 

Proof. We compute below the time evolution of the relative entropy. This is strongly inspired 
from the previous computation, but this takes into consideration the non-conservation of mass 
for the free markers density and the additional dynamics of \i. 

|(m(t)H(i)) = |£ o n(^)(log(^M)-log, + ^)dz 



z>0 



0, (//(/. :)) ( log ( n ^^j ) - log// + ) <1: 



+ / n (t,z)d t log (^-) dz 
Jz>o V rn(t) J 

i 

m(t) 



d z (d z n{t, z) + n(t)n(t, z)) ( log ( ) + vz ) dz 



dt 

where we have used the identity 



z>0 

d . . 
m(t) log v 



'z>0 

We integrate by parts to get 



4(m(t)H(t) + m(t)logi/) = -/ (a g n(t, z) + i*(t)n(t, *)) ( ^"f' 7 + H dz 
dt J z>0 V n (M) / 

-(0,n(i,O)+M(t)n(t,O))log (^r) 

n(t, 2) (0* log n(t, z)) 2 dz + (i/ + fi(t))n(t, 0) 



z>0 



-m(t)/i(t)v - log ( ) 



m{t) J dt 
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We use again the following key identity 



n(t, z) (d z logn(t, z)) 2 dz = I n(t,z) ( d z log n(t, z) + ) dz + 



n(t,0)\ 2 n(t,0) 2 



L>o ' ' ' ' ./, o V '">' ! / " '"f ) 

We end up with the following expression for the dissipation of the corrected entropy, 



£(m(t)H(t)+m(t)logi/) = -/ n(t,z) ( 8 Z log n(t, z) + ^^-) dz 
dt J z>0 V «*(t) / 

n(t,0) 2 



m(t) 



+ {u + (i(t))n(t, 0) - m(t)fi(t)v 



On the first hand, we have that 

" (/ ' ,)! +(i/ + /i(t))n(t,0) -m{t)fi(t)v = (- U<yt, °^ (n(t,0) - m(t)/x) 



m(i) ' \ m(t) 

-m(f) - 

V ' 

- i/) (n(t,0)-m(t)n(t)) 

and on the other hand, we see that 



\ m{t) 



, fn(t,0)\ d , . , \ d 



M (t) 



(n(t, 0) - M (t)) log (^y^) " log (>(*)) ^M*) " (rn(t)) -^m{t) 
(n(i,0)-^))lo"' " (M)) 



— ^ O(t) log //(t) - n(t) + m(t) log m(t) - m(t) - vlog v + M) . 
The last contribution to be reformulated is 

-(ji(t) - u) (n(t, 0) - m(t)fi(t)) = -(p(t) - v) (J^rit) + (1 - m(t))n(t) 

= -(ji(t)-v)(-^/i(t) + (jJi(t)-v)/i(t) 
= ~\^(t) - vf - v(t)(^{t) - vf . 

Combining all these calculations we conclude the proof of Lemma 14.11 □ 

Following the lines of Section \2 . 2. 1 1 we can prove that fi(t) converges to u, the partial mass 
m(t) converges to 1, and the density n(t,-) converges to the stationary state h as t — >• oo. We 
omit the details. 
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5 The higher dimensional case N > 2 

In this section we investigate the possible behaviours of the equation (jl.ip in dimension N > 2 
with the two possible choices (|1.5p and (|1.6p for the advection field. 



5.1 Global existence 

We give the proof of Theorem 11.41 Since many of the arguments are similar to the one- 
dimensional case, we only sketch the proof and focus on the propagation of L p bounds, which 
is the crucial a priori estimate as soon as entropy methods are lacking |23| . 

Let n be a solution of (jl.ip with V • u > and u(t, y, 0) • e z = n(t, y, 0). We see that 



_d 
~dl 



n(t, x) p dx 



H 



-pi Vn{t,x) p ~ l • Vn(t,x)dx 
Jh 

+p Vn(t, x) p_1 • u(i, x) n(t, x) dx . 

Jh 



(5.1) 



On the first hand, we have that 



-pi Vn(t,x) p ~-Vn(t,x)dx 
Ju 



4(p - 1) 



P Jh 



Vn(t,x) p/2 



dx , 



and on the other hand 
P 



P 



- / Vn(t,x) p_1 • u(t,x)n(t,x)dx = - / n(t,x) p (V-u) dx + n(t,y,0) p+1 
l Jh Jh Jx 



dy. 



To estimate the two opposite trends in (|5.ip we use the following Sobolev trace inequality [1] 
and [28j : there exists a constant C r such that for any non-negative / G W ' r we have, 



(7 /(y,0r*dyV A <C r (f \Vf{x)\ r dx\ 



1/r 



(5.2) 



where r* ■ 
estimates: 



^N-^ ■ ^PPlyi n S the previous inequality 



2.s 



/ n(t,y^f s dy<C, 



H 



P J \JH 



Vn(t,x) p/2 



dx 



Vn(t,x) p / 2 n(t,x) s -2 

r*/2 



with / = n s , we obtain the 

r* /r 



H 



n(t, x) 



! — 8 \ 2 



dx 



dx 



(2 — r)r 



We infer that L is the critical space for global existence. Hence we choose 



p\ 2r 
s — - 



2/ 2-r 



N. 



On the other hand, we also require that 



_ r* _ 1 (N - l)r 
~ ~2 ~ 2 N-r 
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A straightforward computation leads to 



2N p + 1 , (2-r)r* 1 

, r s = p + 1 , 



JV + 1 2 ' J 2r N 

Therefore we deduce that 

4- ! n(t, xf dx < - ^ ~ ^ (1 - C\\n(t) II r n) ! Vn(t,x) p / 2 
at J H P Jh 



dx . 



The peculiar choice p = N yields global existence if ||n(0)||£,jv is smaller than some explicit 
threshold as in 1141. 



5.2 Blow-up of solutions in the first case ( 11.51) 

We compute the evolution of the second momentum I(t) = ^ f n \x\ 2 n(t,x) dx as for the 
classical Keller-Segel system (see [30 1 and references therein): 

^1 = NM- [ zn(t,y,0)n(t,y,z)dx. 



H 



Next, define M(t,y) = j z>0 n(t,y,z)dz. Under the assumption d z n(t,x) < for all x G H 
and t > 0, inequality (|2.3ip rewrites 



We deduce that: 



M(t, y) 2 < 2n(t, y, 0) / zn(t, y, z) dz . 

Jz>0 



dK*) / ..m,2 



d/ < NM--\\M(t,y)\\i 2 . (5.3) 
By interpolation there exists a constant C such that 

< Cl{t)^\\M{t,y)\\ 2 ^ . (5.4) 

Indeed we have 



M = [ M(t,y)dy+ [ M{t,y)dy 
J\y\<R J\y\>R 

1 /2 

< CR( N -W([ M(t,yfdy) + iT 2 / \y\ 2 M(t,y)dy 

KJRN- 1 J J R N ~ 1 



< CR^/ 2 \\M(t,y)\\ L2 +R- 2 I(t). 
Optimizing with respect to R we get (|5.4| . Combining (|5.4p and (|5.3p we conclude that the 

JV+l 

solution blows-up in finite time if 1(0) < CM N -^ . 
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5.3 Blow-up in the second case ( 11.61) 

We recall the expression of the advection field in the potential case (|1.6|) : 

(y - v\ z) 



Therefore we have 

dl(t) 
dt 



u(t, x) 

' ' G RiV-l (\y - y'|2 _|_ z 2) 

NM+ I x ■ (n(i, x)u(t, x)) dx 



j^n(t, y', 0) dy' . 



NM 



!)■!! 



y-(y-y') + z 1 



-n(t, y', 0)n(t, y, z) dy dy' dz . 



lz>o (\y - j/f + z 2 ) N/2 
We use a symmetrization trick to evaluate the contribution of interaction: 

[!J '' n(t,y',0)n(t,y,z)dydy'dz = 



y,y> Jz>0 (\y-y>\2 + Z 2f/ 2 



y-y 



I , , >xA frt • [n(t : y',0)n(t,y,z)y-n(t : y,0)n(t,y',z)y') dydy'dz. 

, y > Jz>o (\y-y'\ 2 + z 2 ) 1 

Lemma 5.1. Let f be a smooth positive function. Assume that we have both d z f{x) < and 

Vz > , Vy G M^ -1 , Mh G R N_1 (h • y) (h • d z V y log f(x)) > . (5.5) 

Then for all y, y' G IR^ -1 and for all z > 0, the following inequality holds true: 

(y - v') ■ {f(y', o)/(y, z)y - f(y, o)f(y, z) y ') >\y- y '\ 2 f{ y , z)f(y', z) . 

Proof. Inequality (|5.6p rewrites as follows: 



(5.6) 



(y - y') ■ ( 



f(y,z) 



l 



f(y',z)\ f(y',z) 

y 



l 



f(y,z) 
f(y,o) 



y')>o. 



J(y,o)V f(y',o)J a /(y',o) 

Since d z f(x) < we have both f(y,z) < f(y,0) and f(y',z) < f(y',0) for all y,y',z. Hence 
we are reduced to prove that the vector field 



1 



f(yfi) 



y, 



is monotonic with respect to the y variable. Computing the derivative with respect to y, it is 
straightforward to check that it is monotonic if (|5.5|) is satisfied: 



V, 



f(y,o) 



1 



f(y,o) 



> 




V y f(y,z) V w /(y,0) 



f(y,z) 

Id 



/(y,o) 



<g> y 



V 



s j 9,V„ log /(y, z') dz'^j ® y > , 



in the following matrix sense: A T + ^ > 



□ 
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Under the hypotheses of Theorem 11.51 we assume that conditions (|5.6p - (|5.5p are fulfilled 
for every time of existence. We deduce that 

flI,/) < NM-Uf [ J V ~ y T + *f„ n n(t, y', z)n(t, y, z) dy dy' dz 



dt 2JJy,y>J z > (\y-y>\* + zY /2 

i rr r i 



< NM - - 
2 



y,y' Jz>0 (\y - y'\ 2 + z 2 ) 



iv^nfi) y', z)n(t, y, z) dy dy' dz 



Since \y — y'\ 2 + z 2 < 2\y\ 2 + 2|y'| 2 + z 2 , and n is non-negative, we have 

/{| v |<«M<f,,<*«} 



< NM--IH /? 2 - v ,»(/.//.:)»(/.//.:!d,/d,/<l: 



< ATM — / / n(t,y,z)dy dz 

SB 1 -" ' " 2 



< NM -7—\l / n(i,y,z)dydz , 

where we have used the Cauchy-Schwarz inequality. We have therefore 

<lI(/) < NM-^^(m- [[ n(t,y,z)dydz) 



dt - 4 



{z>^or \y\>f} 



< NM - ^—M 1 + CR~ N ~ 3 I(t) 2 , 
2R _ u.i ^ R 



because R 2 < 9\x\ 2 on > ^ or \y\ > -j}. Optimizing with respect to R, we conclude that 



. JV+l 



the solution blows-up in finite time if 1(0) < CM N ~ 1 , similarly as in Section [5.21 

6 Conclusion 

Here, we have demonstrated that a class of models following [18] exhibit pattern formation 
(either blow-up or convergence towards a non homogeneous steady state) under some condi- 
tions. However we have not answered the main question: do they describe cell polarisation or 
not? Although the one-dimensional case is clear (spontaneous polarisation occurs if the total 
concentration of markers is large enough), the higher-dimensional situation is not so clear. 
Obviously the first model (|1.5p does not exhibit cell polarisation since we can integrate the 
equation (|1.1|) with respect to z, and we obtain for v(t, y) = f z>Q n(t, y, z) dz: 

d t u(t,y) = d yy v(t,y) . 

Thus there is no transversal instability which is the main feature of spontaneous cell po- 
larisation, that leads to symmetry breaking. On the other hand the second model (jl.6p is 
expected to develop symmetry breaking as the tangential component of the advective field 
on the boundary is given by the Hilbert transform of the trace n(t, y, 0) which is known to 
enhance finite time aggregation at least in one dimension of space [12]. However there is no 
clear mathematical distinction between the two models as continuation after the blow-up time 
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appears to be very delicate in a similar context [32], [33j [TTj. It would be very interesting to 
make such a distinction beyond linear analysis as performed in [IB]. We leave it as an open 
question. 
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